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ABSTRACT 

This  thesis  developed  the  ballistics  processor  of  a 
multiple  processor  airborne  tactical  system.   The  multiple 
processor  system  consisted  of  three  INTEL-8080  microcomputers: 
the  executive  processor,  the  navigational  processor  and  the 
ballistics  processor.   The  ballistics  processor  utilized  a 
general  second  order  Runge-Kutta  method  of  integration  of 
the  equations  of  motion  of  unguided  air-to-surface  weapons. 
The  ballistics  processor  computed  sufficiently  accurate  and 
timely  solutions  to  enable  the  executive  processor  to  extra- 
polate an  accurate  release  point  for  the  weapon.   The 
algorithm  permitted  complete  flexibility  in  release  conditions 
and  allowed  a  complete  arsenal  of  air-to-surface  weapons 
currently  carried  on  the  A7-E  aircraft.   The  cost  of  the 
ballistics  processor  using  current  "off  the  shelf"  components 
is  $1635  and  the  entire  tactical  system  was  estimated  to  cost 
$3891. 
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I.   INTRODUCTION 

Military  airborne  tactical  systems  are  used  to  aid  the 
pilot  and  bombardier/navigator  in  their  specific  functions. 
For  an  attack  aircraft,  the  tactical  system  is  used  in  the 
following  manner: 

1.  to  navigate  to  a  given  target 

2.  to  issue  steering  commands  enroute  and/or  to  guide 
the  autopilot 

3.  to  display  to  the  crew  the  current  target  location 
in  relation  to  the  aircraft 

4.  to  solve  ballistics  equations  which  enable  the  crew 
to  release  a  weapon  at  the  appropriate  instant. 

The  presently  operational  systems  employ  a  small 
general  purpose  computer,  such  as  the  IBM  4  PI  system.   In 
this  thesis,  the  ballistics  processor  of  a  multiple  processor 
organization  of  a  typical  tactical  system  for  an  attack  air- 
craft is  presented.   The  system's  architecture  made  use  of 
three  INTEL- 8080  microcomputers,  each  of  which  was  dedicated 
full  time  to  its  assigned  function.   The  purpose  of  the 
multiple  processor  organization  is  to  exploit  the  presently 
emerging  inexpensive  microprocessor  technology  to  solve 
military  tactical  problems.   The  system's  organization  is 
shown  in  Figure  1.   Because  microcomputers  presently  have 
no  hardware  mult iply-and-divide  functions,  the  programmed 
software  mult iply- and-divide  functions  caused  the  microcomputers 
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to  run  on  a  typical  mix  of  operations  approximately  six  times 
slower  than  the  presently  operational  computer  systems.   By 
reorganization  of  the  computations  and  the  use  of  table 
look-up  techniques,  the  described  system  was  capable  of  meet- 
ing the  real-time  requirements  of  an  airborne  weapon  delivery 
system. 

The  ballistics  computer  uses  the  best  estimate  of  the 
position,  dive  angle,  altitude,  and  velocity  of  the  aircraft 
in  its  computation.   It  is  also  provided  with  the  weapon 
identification  (for  the  particular  drag  characteristics) 
and  the  wind  velocities.   From  this  information  and  the 
target's  altitude  (above  Mean  Sea  Level),  the  computer 
determines  the  down  range  travel  and  time  of  fall  of  the 
unguided  weapon  (i.e.,  the  impact  coordinates  in  the  target 
plane) . 

The  accuracy  with  which  a  ballistic  weapon  can  be 
delivered  against  a  target  depends  greatly  upon  the  accuracy 
of  the  sensor- supplied  release  parameters  as  well  as  the 
accuracy  with  which  the  ballistics  equations  are  solved  to 
predict  the  weapon's  down  range  travel. 

Most  aircraft  unguided  air-to-ground  weapons  can  be 
described  as  ballistic  projectiles.   The  only  forces  acting 
on  them  after  release  from  the  aircraft  are  gravity  and 
aerodynamic  drag.   Bullets,  streamlined  bombs,  drogued 
(retarded)  bombs,  cluster  munitions  and  unguided  rockets 
(after  motor  burnout)  are  all  ballistic  projectiles.   Guided 
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weapons  and  weapons  developing  lift  are  not  ballistic  pro- 
jectiles nor  are  they  considered  in  this  development, 

This  thesis  first  describes  the  over-all  tactical  system 
in  order  to  give  the  reader  a  perspective  of  how  the  ballistics 
processor  functions  within  the  entire  weapon  system.   Included 
in  this  background  material  is  a  discussion  of  the  naviga- 
tional computer,  the  ballistics  computer,  the  executive  com- 
puter and  the  intercomputer  communications  network. 

The  basic  calculations  involved  in  solving  the  ballistics 
problem  is  described  next.   The  assumptions  used  in  the  pro- 
blem formulation  are  stated;  the  derivation  of  the  particular 
differential  equations  is  presented;  and  the  application  of 
the  Runge-Kutta  integration  method  is  shown  to  be  suitable 
for  current  operational  airborne  weapon  systems.   The  four 
general  classes  of  weapons  and  how  the  algorithm  handles  each 
integration  is  defined.   The  unique  method  of  compressing  the 
drag  data  for  all  twenty-eight  weapons  is  shown  to  be  a  simple 
curve  fitting  of  the  drag  coefficients  with  respect  to  seme 
standard. 

The  simulation  was  initially  programmed  using  FORTRAN 
for  ease  of  testing  and  debugging.   The  important  phases  in 
programming  were  in  the  DERIV  and  RUNGE  subroutines  where  a 
close  study  of  the  optimum  calculation  methods  was  conducted. 


The  program  was  tested  for  accuracy  of  results  by  comparison 


y 


with  the  published  range  tables  (NAVAIR  01-lC-lT-l). 

Upon  completion  of  the  FORTRAN  implementation,  the  re- 
sultant program  was  translated  into  PL/M,  a  high  level 
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language  compatible  with  the  selected  microcomputer  (INTELLEC 
80) .   The  selected  test  cases  were  verified  correct  in  the 
PL/M  version  of  the  program. 
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I I ,   BACKGROUND 

In  order  to  give  the  reader  a  complete  understanding  of 
how  the  ballistics  processor  functions,  this  section  is 
devoted  to  describing  the  overall  tactical  system  in  which 
the  ballistics  processor  is  an  important  part.   The  first 
section  describes  the  function  of  the  navigational  computer, 
the  sensors  which  provide  data  to  the  computer  and  the 
algorithm  used  to  compute  present  position. 

The  ballistics  computer  will  be  completely  defined  in 
the  remaining  sections  of  this  thesis  and  will  only  be 
briefly  reviewed  here.   Using  the  best  estimate  of  position 
and  velocity  at  a  given  instant,  together  with  information 
of  the  weapon's  drag  characteristics,  the  ballistics 
computer  solves  a  system  of  four  differential  equations  and 
calculates  the  impact  coordinates  of  the  weapon  in  the  target 
plane. 

The  executive  computer  controls  the  display  and  generates 
commands  to  the  autopilot  and  weapon  release  mechanism.   This 
computer  is  interrupted  asynchronously  as  fresh  data  are 
developed  in  the  navigational  and  the  ballistics  computers. 
The  executive  computer  extrapolates  from  the  data  the  appro- 
priate release  instant  of  the  weapon. 

The  intercomputer  communication  net  is  described  next. 
The  star-like  communication  net  provided  a  simple  and  rapid 
means  of  information  transfer. 
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A.   THE  NAVIGATIONAL  COMPUTER 

This  subsystem  is  present  in  all  tactical  systems.   In 
operational  systems,  a  subprogram  (module)  is  periodically 
executed  to  update  the  present  position  by  the  change  in 
position  in  the  last  time  increment. 

In  this  computer,  the  sensor  instruments  provide  digital 
data  to  the  computer  at  specified  time  intervals ,  controlled 
by  a  clock.  There  are  four  sets  of  sensor  instruments  which 
provide  the  data: 

1.  The  Inertial  Navigation  System  (INS) 

2.  The  Doppler  Radar  (DR) 

3.  Air  Data  (AD) 

4.  Radar  Set  (RS) . 

The  INS  system  provides  attitude  (roll  and  elevation)  and 
velocity  increments  in  the  x,  y,  z  directions.   The  DR  is  a 
velocity  sensor  that  uses  the  doppler  principle  for  the 
continuous  measurement  of  aircraft  ground  track  speed  and 
drift  angle.   The  AD  computes  outputs  of  modified  corrected 
static  pressure,  pressure  altitude  and  mach  number  from 
static  and  pitot  tube  pressure  inputs.   The  RS  provides 
target  azimuth,  elevation,  and  range  signals  to  the  naviga- 
tional computer.   A  magnetic  compass   measures  the  aircraft 
heading  in  relation  to  magnetic,  north.   In  order  to  navigate 
accurately  in  this  mode  the  wind  speed  and  direction  are 
determined  by  using  the  "sensed"  velocities  of  the  aircraft. 
Altimeter  measurements  are  also  sampled  in  this  system  to 
update  the  present  position  altitude. 
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There  are  four  basic  modes  in  which  this  system  operates. 
The  crew  has  a  choice  of  using  both  the  INS  (velocities)  and 
DR  outputs  to  minimize  the  Schuler  oscillations.   This  is  the 
most  reliable  and  accurate  mode.   If  the  crew  suspects  that 
either  the  INS  or  the  DR  system  is  malfunctioning,  they  may 
use  only  the  data  from  one  of  the  sources.   The  least 
reliable  and  inaccurate  mode  is  the  last  resort,  when  both 
the  DR  and  the  INS  are  malfunctioning.   The  AD  mode  uses  the 
externally  estimated  wind  vector  and  true  airspeed  of  the  air- 
craft to  determine  the  position  in  latitude  and  longitude. 

In  all  operating  modes,  the  same  basic  technique  of  smooth 
ing  input  data  is  used.   The  least  squares  curve  fitting 
technique  is  used  to  fit  a  quadratic  polynomial  to  the  most 
recent  n  (n  $  16)  data  values.   Use  of  the  Legendre 
orthogonal  polynomial  functions  is  made  in  order  to  make  the 
numerical  calculations  as  accurate  and  rapid  as  possible. 
Smoothing  the  short  time  increment  data  values  will  generate 
less  total  error  in  the  results. 

Using  the  smoothed  data  values,  the  appropriate  calcu- 
lations are  made  to  estimate  the  distance  increments  in  the 
past  time  interval.   As  soon  as  the  new  value  is  generated, 
the  display  of  the  navigational  computer  is  refreshed,  the 
executive  computer  is  interrupted  and  this  most  recent 
estimate  is  transferred  to  the  executive  computer.   If  the 
navigational  computer  is  functional,  the  crew  has  information 
on  the  present  position,  independent  of  the  executive 
computer . 
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The  least  squares  curve  fitting  technique  enables  the 
system  to  diagnose  sensor  malfunction  quite  easily.   If  the 
sensor  generates  fluctuations  in  the  input  data  above  a 
certain  threshold,  then  the  sensor  can  be  diagnosed  as  mal- 
functioning.  Similarly,  if  two  sensors  disagree  above  an 
allowable  threshold,  then  an  automatic  data  analysis  will 
suggest  which  sensor  is  failing.   If  the  crew  concurs  in 
this  diagnosis,  the  operating  mode  of  the  navigational 
system  can  be  changed  to  exclude  the  faulty  sensor. 

B.  THE  BALLISTICS  COMPUTER 

The  ballistics  computer  is  provided  with  the  present 
estimates  of  the  position,  dive  angle,  altitude,  and 
velocity  of  the  aircraft.   It  is  also  provided  with  the  type 
of  weapon  and  wind  conditions.   From  this  information  and 
the  target's  altitude,  the  computer  determines  the  down- 
range  travel  of  the  weapon,  if  released  at  this  instant. 
The  assumption  of  no  guidance  after  release  from  the  air- 
craft is  still  applicable.   However,  the  weapon  may  have  a 
rocket  assist,  or  a  retarded  mode  of  fall. 

C.  THE  EXECUTIVE  COMPUTER 

The  executive  computer  controls  the  main  display  and 
generates  steering  commands  to  the  autopilot  and  firing 
pulses  to  the  weapons  release  mechanisms.   One  of  its 
critical  functions  is  to  extrapolate  from  the  data  generated 
by  the  ballistics  computer  the  weapon  release  point. 
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A  real  time  clock  is  used  to  record  the  time  at  which 
the  asynchronous  data  were  generated  by  the  ballistics 
computer.   The  position  of  the  impact  point  in  the  target 
plane  can  be  described  by  a  quadratic  function  of  time. 
Therefore,  the  real  time  release  point  (adjusted  for  time 
delays)  can  be  extrapolated  very  accurately  from  these  data 
without  additional  extensive  calculations. 

D.   THE  INTERCOMPUTER  COMMUNICATIONS  NET 

In  order  to  keep  the  design  of  the  intercomputer 
communications  net  simple  and  inexpensive,  the  starlike 
network  was  chosen.   The  executive  computer  is  the  center  of 
the  star  and  one  input  port  is  multiplexed  to  the  two  output 
ports  of  the  navigational  and  ballistics  computer  as  shown 
in  Figure  2. 

The  intercommunication  takes  place  in  the  following 
manner:   An  initializing  phase  of  the  program  interrogates 
the  keyboard  to  obtain  data  about  the  initial  position  of 
the  aircraft.   The  latitude  and  longitude  are  entered  in 
degrees  and  minutes  to  the  nearest  1/10  minute;  the 
altitude  is  initialized  at  the  height  above  mean  sea  level. 
These  data  are  made  available  to  the  navigational  computer. 
Also  the  type  of  weapon  and  the  position  of  the  target  (if 
known  in  advance)  is  made  available  to  the  ballistics 
computer.   Similarly,  the  mode  of  navigation  is  determined 
(if  possible  prior  to  flight)  and  relayed  to  the  navigational 
computer. 
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The  execution  phase  of  the  computer  program  makes  use  of 
the  real  time  clock  in  the  executive  computer  which  is 
asynchronously  inrerruping  the  navigational  computer  in 
order  that  sensor  information  be  transferred  via  input  port  2 
to  the  navigational  computer.   When  the  navigational  computer 
has  completed  its  calculations,  it  interrupts  the  executive 
computer  and  transfers  the  calculated  current  position 
sequentially  via  input  port  1  to  the  main  computer.   The 
amount  of  information  sent  is  9  bytes,  and  the  transfer  is 
accomplished  in  about  250  micro-seconds. 

Similarly,  when  the  ballistics  computer  has  completed 
its  calculations,  it  interrupts  the  executive  computer  and 
transfers  its  information  sequentially  via  input  port  1  to 
the  executive  computer.   Only  6  bytes  of  information  are 
sent  in  approximately  170  micro- seconds .   Between  interrup- 
tions, the  executive  computer  updates  the  display  using 
output  port  2.   Calculations  are  continuously  made  to  extra- 
polate the  release  instant  of  the  weapon,  and  the  release 
mechanism  is  activated  when  appropriate  via  output  port  3. 
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III.   BALLISTICS  PROBLEM 

The  equations  of  motion  governing  the  trajectory  of  a 
ballistic  projectile  are  a  simple-appearing  set  of  second- 
order  differential  equations  whenever  the  drag  is  neglected. 
However,  when  a  reasonably  accurate  model  of  the  aerodynamic 
drag  caused  by  the  projectile's  motion  through  the  air  mass 
is  included  in  these  equations  they  are  rendered  nonlinear, 
and  no  general  solution  in  closed  form  has  been  found  that 
is  satisfactory  for  solution  by  an  airborne  fire-control 
computer. 

A.   BASIC  ASSUMPTIONS 

The  mathematical  model  was  chosen  to  represent  reality 
as  accurately  as  possible.   The  armed  forces  publish  range 
tables  for  various  weapons  and  therefore  the  same  model  was 
chosen  as  the  one  used  to  construct  these  tables  (NAVAIR  01- 
1C-1T-1) . 

The  equations  of  motion  were  developed  assuming  the 
projectile  is  a  point  mass  acted  on  only  by  the  force  of 
gravity  and  the  retardation  forces  due  to  air  resistance. 
The  trajectory  was  restricted  to  a  plane  by  ignoring  the 
crosstrail  effects  of  the  wind.   These  crosstrail  effects 
and  the  component  of  wind  which  affects  drag  are  considered 
in  the  executive  computer. 
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The  assumptions  are  summarized  below: 

1.  The  Earth  is  flat  and  nonrotating. 

2.  The  gravitational  attraction  is  constant. 

3.  The  projectile  is  a  point  mass. 

4.  The  projectile  is  not  powered  (except  rockets)  and 
has  a  constant  mass. 

B.   DERIVATION  OF  THE  EQUATIONS 

The  weapon's  fall  is  assumed  to  be  governed  by  a  drag 
force  which  is  proportional  to  the  square  of  the  velocity. 
Newton's  equation  of  motion  becomes 

m  dVt/dt  =  -  g  m  -  ||vt||  2  Cw  Vt/  ||vt|| 

where  the  only  two  forces  acting  on  the  weapon  are  gravity 
and  aerodynamic  drag. 

Vt  -  total  velocity 

m   -  mass  of  the  projectile 

g   -  force  due  to  the  Earth's  gravitation 

Cw  -  drag  coefficient  which  depends  on  the 

characteristics  of  the  weapon  and  may  vary 
with  time  (in  case  the  weapon  has  two  drag 
curves),  and  it  varies  slightly  with  airspeed. 
In  order  to  solve  this  differential  equation,  initial 
conditions,  that  is  the  position  and  velocity,  must  be 
specified.   The  most  convenient  coordinate  system  is  height 
"y"  above  sea  level  and  "x"  downrange  travel,  which  is  zero 
at  the  time  of  release. 
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In  the  x  and  y  component  form,  Newton's  equation  becomes 

dVx/dt  =  -  ||vt||  /m  *  Vx  *  Cw 

dVy/dt  =  -  g  -  ||vt||  /m  *  Vy  *  Cw 

dx/dt   =  Vx 

Dy/dt   =  Vy 

x(0)    =  0  downrange  travel 

y(0)    =  present  altitude 

Vs(0)   =  present  horizontal  velocity 

Vy(0)   =  present  vertical  velocity 

C.   RUNGE-KUTTA  APPLICATION 

The  choice  of  the  "best"  Runge-Kutta  (R-K)  formula  for 
a  particular  problem  was  not  entirely  obvious.   The  lower 
order  formulae  are  simple,  but  require  more  integration 
steps  than  the  more  accurate  higher  order  formulae.   On  the 
other  hand,  the  complexity  of  the  higher  order  formulae 
may  require  so  much  computer  time  that  this  becomes  the 
overriding  consideration.   Generally,  it  was  found  to  be 
safer  to  use  the  lower  order  formulae  when  discontinuities 
in  the  coefficients  of  the  differential  equations  can  be 
expected.   Such  discontinuities  frequently  occur  in  bomb 
trajectories  (e.g.,  drogue  deployment  or  thrust  termination) 
The  second-order  R-K  formula  also  afforded  the  greatest 
flexibility  in  the  distribution  and  size  of  the  integration 
steps.   These  considerations,  along  with  their  simplicity, 
led  to  the  choice  of  the  second-order  R-K  method  as  the  most 
appropriate  for  the  ballistics  algorithm. 
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Let  UI   =  Vx  Fl  =  -1/m  *   |Vt||  *  Vx  *  Cw 

U2   =   y  F2  =  Vy 

U3  =  Vy  F4  =  -1/m  *  ||  Vt  ||  *  Vy  -  Cw  -  g 

The  differential  equations  may  be  expressed  as 

dUi/dt  =  Fi(Ul,U2,U3,U4)         i  =  1,2,3,4 

Therefore,  the  Runge-Kutta  second  order  method  used  in  the 
program  becomes 

mli  =  At  *  Fi(Ul,U2,U3,U4)        i  =  1,2,3,4 

m2i  =  At  *  Fi(Ul  +  cAtmll,U2  +  cAtml2,U3  +  cAtml3, 

U4  +  cAtml4)        c  =  0.7,  i  =  1,2,3,4 
U'i  =  Ui  +  l/2c  *  ((2c  -  1)  *  mli  +  m2i)   i  =  1,2,3,4 

The  procedure  computes  new  values  for  the  variables 
until  the  remaining  time  of  fall,  DTV,  (as  computed  for  a 
vacuum  release) ,  is  less  than  the  value  of  the  next  integra- 
tion time  increment.   The  procedure  then  uses  the  value  of 
DTV  for  the  final  integration  step  size,  and  in  the  calcula- 
tion of  the  down  range  travel  of  the  weapon. 

D.   CLASSES  OF  WEAPONS 

1 •   Streamlined  Bombs  and  Bullets 

The  process  described  above  represents  streamlined 
bombs  quite  well.   After  being  fired,  a  bullet  is  really  just 
a  small,  streamlined  bomb  and  was  handled  as  such.   The  muzzle 
velocity  (3300  feet  per  second)  of  a  particular  gun-bullet 
combination  was  added  to  the  aircraft  velocity. 
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2 .  Drogued  Bombs 

These  bombs  are  more  complex  in  their  trajectory 
modeling.   Their  common  characteristic  is  that  they  are  all 
released  in  a  relatively  low-drag  configuration.   At  some 
time  after  release,  they  deploy  vanes  or  a  parachute  which 
greatly  increases  the  total  drag.   To  compensate  for  this 
change  of  drag  characteristics,  the  initial  integration  step 
size  was  chosen  to  coincide  with  the  deployment  of  the  high 
drag  devices.   Thereafter,  there  was  no  difficulty  in  calcu- 
lating the  trajectory  of  the  drogued  weapons. 

3.  Cluster  Bombs 

These  weapons  are  released  as  large,  low-drag 
containers.   At  a  predetermined  point  in  their  trajectory, 
which  is  time  dependent,  they  use  some  method  to  dispense 
smaller  weapons.   These  smaller  weapons  are  of  higher  drag 
than  the  container.   The  ballistics  computer  will  predict 
the  impact  point  of  the  center  of  the  pattern  rather  than 
individual  positions  of  the  bomblets.   The  method  of  applying 
the  algorithm  is  to  calculate  the  trajectory  of  the  container 
to  the  point  where  it  dispenses  the  small  weapons;  switch 
to  a  different  drag  curve  which  describes  the  motion  of  the 
pattern  center;  and  integrate  down  to  the  target  altitude. 
As  in  the  retarded  bombs,  the  initial  integration  step  was 
chosen  to  coincide  with  the  act  of  dispensing  the  small 
weapons . 

^ ■   Unguided  Rockets 

Rockets  present  several  unique  complexities.   They 
change  weight  during  flight;  they  have  thrust  as  well  as 
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drag;  and  they  slew  around  after  firing,  because  their 
launches  may  not  line  up  with  the  aircraft  direction  of 
flight.   The  algorithm  treats  unguided  rockets  like  a 
cluster  bomb  with  the  first  stage  having  an  average  thrust 
as  well  as  an  average  thrust  time.   This  provides  sufficient 
accuracy  and  saves  computer  storage  that  would  be  taken 
storing  the  thrust  profile  of  each  rocket. 

E.   COMPRESSION  OF  DRAG  DATA 

The  behavior  of  a  particular  weapon's  drag  coefficient, 
Cd,  as  a  function  of  Mach  number  must  be  known  to  calculate 
impact  ranges  for  that  weapon.   Normally,  this  information 
is  given  in  tabular  form  with  some  weapons  having  more  than 
100  points  in  the  table.   Most  aircraft  have  the  capability 
to  carry  many  different  types  of  weapons.   It  is  desirable 
to  store  all  the  Cd  information  of  an  aircraft's  weapon 
repertoire  in  the  airborne  computer  to  avoid  loading  the  Cd 
data  for  each  different  combination  of  weapons.   This 
requires  storing  drag  coefficient  data  in  a  more  efficient 
form  than  a  table. 

The  following  factors  influence  the  choice  of  a  scheme 
to  approximate  the  drag  coefficients  for  a  given  set  of 
weapons : 

1.  Computer  storage 

2.  Computer  computation  time 

3.  Delivery  envelope  of  each  weapon 

4.  Allowable  downrange  and  time-of-fall  errors. 
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It  would  be  ideal  to  use  one  general,  fairly  powerful 
expression  to  fit  all  the  drag  tables  for  the  weapon 
repertoire  of  a  given  aircraft.   This  is  quite  difficult  to 
accomplish  because  weapons  like  guns  and  rockets  require 
Cd  values  for  Mach  numbers  much  higher  than  free-fall  bombs. 
As  a  consequence  of  this  generally  acceptable  tabling  pro- 
cedure, extra  storage  is  spent  on  weapons  that  do  not  require 
the  extra  capability. 

A  significant  saving  can  be  realized  when  it  is 
recognized  that  several  free-fall  weapons  have  drag 
coefficient  curves  that  differ  only  by  a  multiplicative 
factor.   The  Mk  80  low-drag  series  are  one  such  group  of 
weapons  with  the  Mk  84  drag  curve  conventially  taken  as  the 
reference.   Similarly,  the  Mk  106,  CBU  bomblet,  and  Sadeye 
bomblet  are  another  group  with  the  Garve  drag  curve  as  the 
reference.   The  specific  region  of  interest  of  many  more 
drag  curves  can  be  approximated  by  multiplying  these  two 
reference  drag  curves  by  a  factor  and  then  translating  them 
along  the  Cd  and  M  axes  until  they  match  the  original  drag 
curves . 

To  further  aid  in  the  explanation  of  approximating 
weapon  drag  curves  by  this  method,  let  Cd(M)  and  M  be  the 
drag  coefficient  and  Mach  number  of  a  weapon  drag  curve. 
Let  Cd (M)  be  the  drag  coefficient  of  a  reference  drag  curve 
and  let  A,  B  and  DM  be  the  multiplicative  factor,  drag 
coefficient  translation,  and  Mach  number  translation, respec- 
tively, that  when  applied  to  the  reference  drag  curve  will 
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approximate  the  particular  weapon  drag  curve  of  interest. 
That  is: 

Cd(M)  *  A  *  Cd(M  +  DM)  +  B 

A  hypothetical  weapon  drag  curve,  Cd,  and  reference  drag 
curve,  Cd,  are  shown  in  Figure  3.   For  simplicity  of  discus- 
sion, the  Mach  number  intervals  (M2  -Ml)  and  (M4  -  M3)  are 
equal.   Suppose  the  Cd  curve  is  multiplied  by  an  A  so  that 
the  Cd  and  the  A  *  Cd  curves  are  nearly  identical.   For  the 
Cd  and  A  *  Cd  curves  to  be  a  congruent,  the  A  *  C d  curve  must 
be  translated  along  the  Mach  axis  an  amount  equal  to  DM,  and 
along  the  Drag  Coefficient  axis  an  amount  equal  to  B. 

For  many  weapons,  the  above  stated  equation  will  not  give 
an  accurate  fit  over  a  wide  interval,  of  Mach  numbers;  instead, 
A,  B  and  DM  should  be  chosen  so  that  the  Cd  fit  is  best  for 
the  Mach  number  interval  that  corresponds  to  the  delivery  en- 
velope for  a  particular  aircraft. 

Considering  the  weapons  currently  implemented  on  a  typical 
attack  aircraft,  (A7-E),  both  the  Mk  84  and  the  Garve  reference 
drag  curves  are  necessary  to  handle  the  free-fall  weapons. 
Instead  of  extending  these  curves  to  handle  guns  and  rockets 
that  have  a  much  larger  Mach  number  range,  a  rocket-gun 
reference  drag  curve  is  used. 

F.   CURVE  FITTING  THE  DRAG  COEFFICIENTS 

In  the  previous  section,  the  reference  curves  needed  and 
how  they  were  used  to  approximate  weapon  drag  curves  for  use 
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in  the  algorithm  was  presented.   In  this  section  the  reference 
drag  curves  versus  Mach  number  is  represented.   The  form  of 
the  function  used  will  minimize  computer  storage,  logic,  and 
computation  time  for  the  allowable  error  in  downrange  and 
time  of  fall  for  the  different  weapons.   It  was  decided  that 
each  of  the  three  reference  curves  be  divided  into  three 
regions  with  Cd  expressed  as  a  second-order  polynomial  in 
Mach  number  for  each  region. 

Cd(M)  =  aO  +  al  *  M  +  a2  *  M2 

The  coefficients,  aO,  al,  and  a2  depend  on  the  region  the 
Mach  number  M  is  in.   The  drag  coefficient  for  a  particular 
weapon  is  expressed  as: 

Cd(M)  =  A  *  (aO  +  al  *  (M  +  DM)  +  a2  *  (M  +  DM)2)  +  B 

As  a  further  step  to  save  computer  storage  and  reduce  the 
computation  time  for  the  algorithm,  another  compression  of 
drag  data  was  undertaken.   The  computation  C  *  (PI/8  *  Cd) , 

CKDG,  is  made  in  the  Runge-Kutta  portion  of  the  algorithm 

2 
where  C  is  the  bomb  factor,  d  /W,  where  d  is  the  weapon 

diameter  in  feet  and  W  is  the  mass  in  slugs.   One  constant 

per  weapon  in  computer  storage  and  two  multiplications  per 

integration  step  is  saved  by  expressing  CKDG  as 

CKDG  =  CF  *  CKDG  +  DKG 

where    CF  =  A  *  (d/d) 2  *  W/W 


31 


CKDG  =  d  /W  *  PI/8  *  Cd 

DKG  =  d2/W  *  PI/8  *  B 

In  the  above  equations,  d  and  W  are  the  assumed  diameter 
and  mass  for  the  reference  curve  Cd.   The  final  form  of  CKDG 
in  the  algorithm  is 

CKDG(M)  =  CF  *  (bO  +  bl  *  (M  +  DM)  +  b2  *  (M  +  DM)2)  +  DKG 

Table  1  gives  a  list  of  bO,  bl,  and  b2  for  each  region  of  the 
three  reference  curves  and  the  values  of  CF,  DM,  and  DKG,  are 
tabulated  in  Table  2. 

Thrust  appears  in  the  algorithm  in  the  total  drag 
function 

HH  =  TH/V  -  RHO  *  CKDG  *  V 

where   TH  =  thrust/mass 

RHO  -  air  density 

V  =  weapon  velocity 

and  CKDG  was  previously  defined  in  the  above  equations. 
Figure  4  shows  thrust  versus  time  for  a  typical  rocket 
approximated  by  a  constant,  thrust  =  total  impulse/Tl, 
where  Tl  is  decreased  from  actual  motor  burn  time  in  order 
to  make  thrust  equal  to  the  average  thrust  of  region  2.   The 
mass  of  the  rocket,  W,  obviously  decreases  during  rocket 
motor  burn  time.   However,  it  is  adequate  for  W  to  be 
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TABLE  1.  Reference  Drag  Curve  Coefficients  and  Cuts. 
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approximated  by  a  constant,  by  varying  W  until  dov/nrange 
error  is  minimized.   The  second  stage  of  the  rocket,  for  time 
after  weapon  release  greater  than  Tl,  is  free  fall  and  TH 
will  be  zero. 
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IV.   FORTRAN  SIMULATION 

A.   BOMB  DROP  SIMULATION 

The  initial  ballistic,  trajectory  algorithm  for  digital 
airborne  fire  control  computers  used  in  the  simulation  of 
Traj ,  Deriv,  and  the  Runge  routines,  was  developed  by  the 
Boeing  Company,  Seattle,  Washington.   The  research  was 
completed  under  a  Naval  Weapons  Center  contract  in  September 
1970. 

As  previously  mentioned,  the  Second-Order  Runge-Kutta 
integration  scheme  was  chosen  because  of  its  flexibility 
and  simplicity.   The  results,  as  evidenced  in  Tables  3,  4 
and  5,  compared  favorably  to  the  published  data  (NAVAIR 
01-lC-lT-l),  with  mean  errors  of  0.101  for  the  down  range 
travel.   The  low-drag  bombs  showed  relatively  little 
sensitivity  to  the  step  size,  so  long  as  the  maximum 
integration  step,  (DMAX) ,  was  less  than  or  equal  to  5 
seconds  and  the  normal  Initial  Integration  Step,  (DTI), 
was  less  than  or  equal  to  3  seconds.   However,  the  dual 
stage  weapons  (IDN0>17)  did  show  a  large  sensitivity  to 
the  step  size  used  during  the  integration  scheme.   The 
accuracy  desired,  required  the  DTI  value  to  be  consistently 
less  than  or  equal  to  2  seconds  with  the  average  value  of 
1.05  seconds,  for  the  final  eleven  weapons. 

The  ultimate  verification  of  the  step  size  for  the 
Runge-Kutta  (Deriv  also)  method  was  conducted  by  varying 
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TABLE  3.  Sample  Errors  in  Impact  Range  Calculations  for  the 
General-Purpose  Bomb  Mk  82 /Snakeye /Unre tarded . 


Release 

Conditions 

ImDac t  Ranee,  ft 

Impac  t  Rang  e 

Velocity,  - 

altitude  , 

Angle  , 

Ballistic 

Algor i  thm 

error , 

kno  t  s 

ft 

deg 

Tables 

ft 

400 

2500 

0 

8039 

8039 

0 

5000 

-10 

9158 

9159 

+  1 

12000 

-30 

10536 

10537 

+  1 

12000 

-45 

7530 

7531 

+  1 

450 

2500 

0 

8998 

899  5 

-3 

5000 

-10 

10025 

10023 

-2 

12000 

-30 

11340 

11341 

+  1 

12000 

-45 

7998 

7999 

+  1 

500 

2500 

0 

9946 

9940 

-6 

5000 

-10 

10843 

10839 

-4 

12000 

-30 

12067 

12067 

0 

12000 

-45 

8406 

8407 

+  1 

TABLE  4.  Sample  Errors  in  Impact  Range  Calculations  for  the 
Special-Purpose  Bomb  Mk  40/Mcd  0/Conical  Tail. 


Release 

Conditioi 

is 

Imnact  Ranee,  ft 

Impact  Range 

Velocity  , 

Altitude  , 

Angle  , 

Ballistic 

A.1  gorit  hm 

error  , 

knots 

ft 

de  2; 

Tables 

ft 

400 

2500 

0 

8305 

8301 

-4 

5000 

-10 

9434 

9432 

-2 

12000 

-30 

10847 

10846 

_  i 

12000 

-45 

7691 

7690 

-  1 

450 

2500 

0 

9329 

9324 

-5 

5000 

-10 

10342 

10339 

-3 

12000 

-30 

11670 

11669 

-1 

12000 

-45 

8160 

8159 

-1 

500 

2500 

0 

10348 

10342 

-6 

5000 

-10 

11200 

11197 

-3 

12000 

-30 

12401 

12402 

+  1 

12000 

-45 

8559 

8559 

0 
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TABLE  5.  Sample  Errors  in  Impact  Range  Calculations  for 
Practise-Purpose  Bomb  Mk  76/  with  Lug. 


the 


Release  Conditions 

Impact  Ranee,  ft 

Imp  act  Range 

Velocity  , 

Vlt i tude  , 

Angle  , 

Ballistic 

Algo  r i  thm 

error 

kno  t  s 

ft 

deg 

Tables 

ft 

400 

2500 

0 

7633 

7617 

-16 

5000 

-10 

8719 

8710 

-9 

12000 

-30 

9991 

9992 

+  1 

12000 

-45 

7232 

7235 

+  3 

450 

2500 

0 

8463 

8459 

-4 

5000 

-10 

9493 

9487 

-6 

12000 

-30 

10712 

10720 

+  8 

12000 

-45 

7669 

7674 

+  5 

500 

2500 

0 

9268 

9255 

-13 

5000 

-10 

10195 

10192 

-3 

12000 

-30 

11337 

11340 

+  3 

12000 

-45 

8038 

8039 

+  1 
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the  time  interval  from  one  to  five  seconds.   The  TRAJ 
integration  scheme  was  simultaneously  tested  using  several 
release  conditions  to  demonstrate  the  model's  flexibility  in 
simulating  bomb  loft,  level  flight  release  and  dive  bombing 
conditions.   The  emphasis,  at  this  point,  was  on  the  number 
of  functional  calls  to  RUNGE  for  the  projectile  to  reach  the 
target  altitude.   The  results,  by  altering  At,  showed  changes 
ranging  from  0.00°6  to  0.44%  for  the  time  of  fall,  and  0.00% 
to  0.46%  for  the  downrange  travel,  with  a  substantial  decrease 
in  the  number  of  calls  to  RUNGE  also  evidenced  (Table  6) . 

B.   OVERALL  SIMULATION 

The  PRINT  routine  was  a  means  of  simulating  the  various 
cockpit  displays  in  the  tactical  attack  aircraft.   The  output 
consisted  of  the  release  altitude,  dive  angle,  and  airspeed, 
as  well  as  the  resulting  downrange  travel  and  the  time  of 
fall  of  the  weapon. 

The  INPUT  routine  simulated  the  collecting  and  passing 
of  the  sensor  data  to  the  various  processors.   The  sensors 
included  in  the  attack  aircraft  are  the  Inertial  Navigation 
System,  Doppler  Radar,  Air  Data  and  the  Radar  Set.   Specifi- 
cally, the  velocity,  dive  angle,  altitude  and  the  identifi- 
cation number  of  the  weapon  (usually  provided  via  the  armament 
panel)  are  required  for  the  execution  of  the  program. 

The  SETDAT  routine  completed  the  interface  between  the 
executive  computer  and  the  ballistics  computer.   Its  purpose 
was  that  of  assigning  the  constants  their  specific  nonvarying 
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values  and  reinitializing  the  variables  to  zero  at  the  start 
of  each  new  weapon's  integration.   The  routine  also  performed 
various  calculations  including  the  determination  of  the 
velocity  in  the  x  and  y  directions,  the  conversion  of  raw 
velocity  to  feet  per  second  and  the  calculation  of  the  release 
angle  in  radians. 

The  weapon  dependent  constants  required  for  the  algorithm 
were  stored  in  DECODE.   Also,  DECODE  contains  the  necessary 
logic  to  load  the  correct  constants  for  the  particular 
weapon  the  crew  has  selected.   After  the  first  pass,  there 
is  no  need  to  go  through  DECODE  before  entering  the  algorithm 
for  a  trajectory  calculation.   The  actual  scheme  used  for 
DECODE  in  reality  will  depend  on  the  particular  aircraft  and 
its  fire  control  requirements  and  how  the  algorithm  is  inter- 
faced in  the  airborne  computer. 
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V.   PL/M  IMPLEMENTATION 

To  implement  the  model  on  the  INTEL  8080  microcomputer, 
the  floating  point  arithmetic  package  was  used  to  maintain  a 
four  hexadecimal  digit  accuracy  during  calculations.   The 
floating  point  variable  therefore  occupies  three  eight-bit 
words:  two  words  for  the  variable  mantissa  and  one  word  for 
the  mantissa  sign  bit  and  the  variable  exponent.   The  man- 
tissa of  the  variable  is  viewed  as  a  normalized  binary  number 
with  the  left  most  digit  always  one. 

In  order  to  manipulate  these  variables,  an  arithmetic 
subroutine  library  was  supplied  by  Professor  Kodres.   At 
present,  there  exist  routines  that  add,  subtract,  multiply, 
divide  and  take  the  square  roots  of  these  floating  point 
variables . 

The  actual  PL/M  version  of  the  ballistics  program  was 
translated  directly  from  the  working  FORTRAN  version.   All 
the  associated  flow  charts,  process  graphs,  program  listings, 
and  the  simulated  output  can  be  found  in  the  appendix. 

The  PL/M  routines  have  been  compiled  and  tested.   The 
results  compare  very  favorably  with  the  FORTRAN  results. 
The  PL/M  verification  of  the  test  parameters  is  included  in 
Table  7. 
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TABLE  7.  PL/M  Test  Verification  (time  interval,  MK  83) 


RELEASE 
CONDITIONS 

BALLISTIC 
TABLES 

— 
DT  =  3.0 

COMPUTAT3 

ON  TIME 

Trial  1 

Trial  2 

10° 
350  kts 
1000  ft 

11.71 

6701, 

11.69  -.171 

(4) 
6694.  5-.  093 

1.2 
(sec) 

1.1 
(sec) 

10° 

350  kts 
3000  ft 

17.29 
9826. 

17.28  -.053 
(6) 

9  8  2  0  .  5  -  .  0  5  3 

1.8 

1.7 

10° 
650  kts 
500  ft 

13.89 
14192. 

13.82  -.49% 

(5) 
14140.  -.363 

1.  3 

1.2 

10° 
650  kts 

3000  ft 

20.85 
2  0  8  5  9. 

20.81  -.193 
(7) 

2  0  8  2  0  .  -  .  1 8  3 

2.  2 

2.1 

0° 
300  kts 
3500  ft 

14.85 

7  3  8  2. 

14.83  -.133 

(5) 
7  3  79.  -.033 

1.6 

1.4 

0° 
300  kts 
15000  ft 

30.98 

1516  4. 

30.95  -.09  3 

(11) 

1 5 1 5  6  .  -  .  0  4  3 

3.2 

3.0 

0° 
450  kts 
5000  ft 

17.82 
1312  7. 

17.81  -.053 
(6) 

1 5 1 2  2  .  -  .  0  3  3 

1.  8 

1.7 

-20° 
400  kts 
5500  ft 

12.81 
7  9  6  3. 

12.80  -.0  73 

(5) 
7963.  0.003 

1.2 

1.1 

-20° 
550  kts 
7000  ft 

13.50 
113  8  7. 

13.47  -.223 

(S) 
11388. +. 013 

1.  2 

i.  2 

-40° 
450  kts 
5000  ft 

8.18 
46  89. 

8.18  0.003 
(5) 
4  6S9.  0.003 

0.9 

1.0 

-60° 
400  kts 
70  0  0  ft 

9.61 
5194. 

9.60  -.103 
(4) 
319  5.  +.023 

0.3 

0.  8 
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VI.   CONCLUSIONS 

The  objective  of  this  thesis  was  to  demonstrate  a  micro- 
computer based  digital  system  which  will  perform  the  same 
functions  as  the  attack  portion  of  a  presently  implemented 
airborne  tactical  system.   The  main  purpose  of  the  project 
was  to  verify  the  feasibility  of  implementing  the  most  com- 
plex arithmetic   task  of  the  airborne  tactical  computer  by 
use  of  an  INTEL  8080  microcomputer  at  an  estimated  cost  of 
$3891  for  the  entire  tactical  system. 

Because  of  the  rapidly  decreasing  hardware  costs  of  the 
microcomputer  and  significant  increases  in  processing  speeds 
which  the  microcomputer  industry  is  presently  realizing, 
these  new,  innovative  systems,  which  are  in  all  ways  equivalent 
to  the  presently  available  avionics  systems,  can  be  developed 
at  a  hardware  cost  which  is  nearly  10  times  smaller.   The 
software  development  costs  are  also  likely  to  be  smaller 
because  the  system  can  be  subdivided  into  smaller,  independ- 
ently operating  subsystems,  each  of  which  will  perform  a  well 
defined  task. 

The  lack  of  hardware  multiply-and-divide  operations  in 
the  microcomputer  is  compensated  by  careful  process  analysis 
to  reduce  the  number  of  multiply-and-divide  operations.   This 
allows  three  ,:slow"  (in  the  computer  sense)  computers  to 
solve  the  same  problem  that  one  "fast"  computer  is  presently 
solving,  at  a  great  hardware  cost  reduction.   The  necessity 
for  careful  analysis  of  the  problem  may  even  disappear  soon, 
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because  the  microcomputer  manufacturers  are  likely  to 
incorporate  the  hardware  multiply-and-divide  operations. 
This  thesis  demonstrates  the  feasibility  of  using  the  micro 
computer  as  a  systems'  building  block  for  airborne  tactical 
systems . 


APPENDIX  A   PROGRAM  FLOW  CHARTS  AND  PROCESS  GRAPHS 


Symbols  Used  In  Flow  Charts 


(        ) 


Start/Stop 


Output 


Card  Input 


Processing  Symbol 
(Arithmetic   Operation) 


Conditional  Symbol 


Subroutine  Call 


Return  Statement 


O 


Program  Entry  Point  (label) 
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Symbols  Used  in  Process  Graphs 

input 


output 


1.   Bite  Operations  (8  bits): 


ADD 
MULTIPLY 

SUBTRACT 

DIVIDE 

ASSIGNMENT 


2.   Floating  Point  Operations: 


ADD 


MULT 


DIV 


SQRT 


ADD  (Subtract)  Routine 


MULTIPLY  Routine 


DIVIDE  Routine 


SQUARE  ROOT  Routine 


3.   Subroutine  Calls 


)ERIV 


First  Derivative  Routine 
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MAIN    PROGRAM 


c 


START 


I 


J 


SET=    -1 


I 


SETDAT 


5 


SWITCH=    1 


PRINT 


5 


O 


INPUT 


c 


STOP 


J^ 


(T) 


IIDN0=1 


KIDNO 


STEP    1 


<y 


I DNO=NUMI D ( I I DNO ) 
SET    =    0.0 


51 


MAIN    PROGRAM    (PART    2) 


SETDAT 


5 


DECODE 


1 


IANG=1 

STEP  1 

KDEG 

<h 


SL 


THETA=DEG(IANG) *RAD 


IY=1 

STEP  1 

kiy 

<h 


2 


SET   =1.0 


5 


<] — (4 


SETDAT 
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MAIN    PROGRAM    (PART    3) 


SWITCH=    -1. 


PRINT 


TRAJ 


£>(5 
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INPUT    SUBROUTINE 


C 


ENTRY 


} 


v 


VELOCITY 


2 


ALTITUDE 


v 


DIVE  ANGLE 


v 


WEAPON  ID 
NUMBER 


54 


PRINT    SUBROUTINE 


C 


ENTRY 


I 


3 


GLOBAL 
VARIABLES 


OUTPUT 
DATA 


2. 


HEADING 
INFORMATION 


55 


SETDAT    SUBROUTINE 


c 


ENTRY 


I 


) 


GLOBAL 
VARIABLES 


7 


G=    32 

.174 

RAD= 

.0174533 

**A=    0.7 

AA=    0 

5/A 

YT=    0 

0 

VYK=    - 

-5.0 

FRACT= 

=    0.5 

RETURN 


CFORMl=  0.0 
CFORM2=  0.0 
DM1=  0.0 
DM2=  0.0 
DKG1=  0.0 
DKG2=  0.0 
VMUZ=  0.0 
VE=  0.0 
SL=  0.0 
DMAX=  5.0 
***ITYPE=  -1 
FN=  0.0 
IBOTH=  1 
**T=  0.0 


5 


X±J 


1 


U=    VKTS(IV) *1.6878 

DEL=    ATAN(VE,U) 

V=SQRT (U*U+VE*VE) 

VXA=    (V+VMUZ) *COS(THETA-DEL) 

VYZ=    (V+VMUZ) *SIN(THETA- DSL) 


i 


RETURN 


RETURN 


*  SET  VALUES  FOR  PL/M  ARE  1,  2,     3 
**  T  IS  TM  AND  A  IS  A   IN  PL/M. 
***  ITiTE  HAS  A  VALUE  3  FOR  -1  IN  PL/M. 
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DECODE    SUBROUTINE 


(  ENTRY  ) 


GLOBAL 
VARIABLES 


c 


4 


IDNO 


) 


0 — &E 


G> 
G> 


REF=    4 
DKG1=    2.55E-03 
DTI=    3.0 


£> 


IREF=    4 

DKG1=    6.29E-03 

DTI=    3. 


G> — t>P 


REF=    4 
DKG1=    4.01E-03 
DTI=    3. 


-£>jlREF=    2 
DMAX=    3 . 

CF0RM1=. 3.92E-03 
DKG1=    2.75E-03 
DTI=    2. 


-IX 


IREF=    2 
DMAX=    3. 

CF0RM1=    3. 91 E- 03 
DKG1=    6.36E-03 
DTI=    3. 


31 


31 


31 


£>n 


IREF=    4 

DMAX=    2 

DKG1=    0. 

021266 

DTI=    1. 

IREF=    1 

CFORMl- 

2.5704 

DTI=    3  . 

o© 


-o 


3] 
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DECODE    SUBROUTINE     (PART    2) 


G^H 


12 


14 


IREF=    4 
DMAX   =    3. 
DKG1-    9.77E-03 
DTI-    2. 


■O 


IREF=    1 
CFORMl-    2.064 
DTI-    3. 


GG^ 


I>REF  =    1 

CFORMl    =    1. 34  31 

DTI=    1. 


tx 


IREF-    1 
CF0RM1=    1.21 
DTI-    3. 


© ^ 


IREF=    1 
CFORM1    =1.0 
DTI-    3. 


€> 


IREF-  1 
CFORMl^  3.12 
DKG1=  -1.22E-03 
DTI-  3. 


t>( 


31 


31 


31 


31 


31 


15 


O 


IREF-    1 
DMAX-    3. 
CFORMl-    3.4972 
DTI-    2. 


£> 


IREF-    ] 
CFORMl-    1.605 
DTI-    3- 


17 


~t> 


IREF-    4 
DMAX-    3. 
DKG1-    7.3  3F-0  3 
DTI-    1 . 


31 
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DECODE  SUBROUTINE  (PART  3) 

18 


19 


20 


21 


O 


IREF=    1 

ITYPE=    1 

IBOTH=    2 

DKG1  =    7.33E-03 

CFORM2=    1.69E-02 

DM2=    0.38 

DKG2  =    0.17166 

DS=    0.6617 

SL=    -0.000269 

DTI=    2. 


t> 


1> 


t> 


€>0 


IREF=    1 
ITYPE=    1 
I BOTH =    2 
CFORMl=    2.0754 
CFORM2=    0.2217 
DS=    4.267 
DTI=    1. 5 


IREF=  1 
I TYPE-  1 
IBOTH=  2 
CFORMl=  2.2973 
CM1=  0.32 
DKG1  =  8.17E-03 
CFORM2=  1.11E-02 
DM2=0.41 
DKG2=0. 16885 
DS=  4.06 
DTI=  2.0 


31 


IREF=  1 
I TYPE =  1 

r 

I 

IBOTH=  2 

CFORMl=  2 

2404 

CFORM2=  0 

1178 

DS=  4.0 

DTI=  1.62 

31 
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DECODE    SUBROUTINE     (PART    4) 

22 


0- 


0 


24 


25 


26 


27 


28 


■o 


IREF=  1 
ITYPE=  1 
IBOTH=  2 
DKG1=  9.77E-03 
CF0RM2=  2. 3 IE- 02 
DM2=  0.38 
DKG2=  0.23287 
DS=  0.679 
SL=  -0.000303 
DTI  =  1.622 


■Ol 


CH 


IREF=  3 
ITYPE=  2 
CFORM1-  0.82 
CFORM2=  1 . 0 
FN=  1746. 
DS=  1.4225 
DTI=  1. 


•O 


IREF=  4 
ITYPE=  0 
DS=  0.98 
DKG2=  1.48 
DTI=  0.31 


C> 


IREF-  4 
ITYPE=  0 
DS=  0.89 
DKG2=  2. 
DTI=  0.22 


£> 


IREF=  4 
ITYPE=  0 
DS  =  0.89 
DKG2=  2.70 


-Ol 


IREF=    2 
ITYPE=    2 
CFORM1-    0.1514 
CFORM2=    0.1514 
DS  =    0.5 
DTI=    0.8 


31 


IREF=  3 

r 

1 

DMAX=  1.5 

CF0RM1=  2.9964 

DKG1  =  -0.014992 

VMUZ=  3300. 

DTI=  0.5 

o© 


31 


31 


31 


31 


31 
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DECODE    SUBROUTINE    (PART    5) 


©— -a 

Q  2JO 1 ^Tref         V  i 


& 


33 


34 


0 


■O 


CC(1,  1,1)  = 
CC(1,  2,1): 

CC ( 1 , 3 , 1 ) = 
CC  ( 2  , 1 , 1 )  = 
CC(2,2,1)= 
CC(2,3,1)= 
CC(3,1,1)= 
CC ( 3 , 2 , 1 ) = 
CC ( 3 , 3 , 1 ) = 
CT(1,1)=0. 
CT(2,1)=Q. 


:1.6E-03 

:0.0 

0.0 

4.7E-02 

-0.1097 

6.7E-02 

-0.1164 

0.2176 

-9.7E-02 

834 

977 


i> 


CC  ( 1 , 1 , 1 )  = 

.1041 

CC  (1,2,1  )=■ 

-.2303 

CC(1, 3,1)= 

.1676 

CC(2,1,1)=- 

CC(2,2,1)= 
CC(2,3,1)=- 

-.1940 

.4015 

-.1646 

CC ( 3 , 1 , 1 ) = 

7.3E-02 

CC(3,2,1)= 
CC(3,3,1)= 
CT(1,I)=  1 
CT(2,1)=  1 

-2.0E-02 

2.44E-03 

032 

30 

-oT  ret 


URN 


O© 


CCU,  1,IB0TH)  = 

3.5350 

CC(1,2,1B0TH)= 

-3.3478 

CC(1,3,IB0TH)= 

2.8726 

CC(2,1,IB0TH)= 

11.2617 

CC(2,2,IBOTH)-- 

-27.4162 

CC(2,3,IBOTH)= 

21.7308 

CC(3,1,IB0TH)= 

-23.7915 

CC(3,2,IBOTH)= 

44.2607 

CC(3,3,IBOTH)= 

-14.4996 

CT(1,IB0TH)=  0 

622 

CT(2,IB0TH)=  0 

885 

f 

I 

o© 


61 


TRAJ    SUBROUTINE 


i 


D=    DS+SL*U 


£ 


C 


ENTRY 


i 


3 


CF= 

CFORMl 

DM= 

DM1 

DKG 

=  DKG1 

MSTG=  1 

X= 

0.0 

T= 

0.0 

VX= 

VXA 

VY= 

VYA 

TH  = 

FN 

YA= 

Y 

Y= 

ALT(IY) 

5 


RUNGS 


5 


DTV=    1/G*(VY+SQRT(VY**2+2.*G*V) ) 
D=    DTI 


(F) 


MSTG=    2     I 


TEH     ' 

MSTG=6   TABLE=27 


o> 


*    ITYPE=    -1    IN    FORTRAN 
3    IN    PLM 


-O 


_i 


D=    DMAX 


MSTG=    1 


M>© 


PLM 
MSTG=  0 
TABLE 1=  0 
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TRAJ    SUBROUTINE    (FART    2) 


©- 


(7>o 


DKG=    DKG2 
DM=    DM2 
CF=    CFORM2 
TH=    0.0 


o 


■o 


DKG=    DKG2 
DM=    DM2 
CF=    CF0RM2 
TH  =    0.0 


D=    DTV 


I  PLM  TABLEl=  0 
MSTG=  1  !  MSTG=  0 


1 


RUNGE 


DTV=1/G* (VY+SQRT(VY**2+2.*G*Y) ) 
T=    T+DTV 
X=    X+DTV*VX 
SWITCH=    0.0 


I 


PRINT 


RETURN 
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RUNGE    SUBROUTINE 


C 


ENTRY 


I 


3 


GLOBAL 
VARIANCES 


I 


AD=    A*D 
YO=    Y 
VXO=    VX 
VYO=   VY 


2 


RHO 


=2.4E-04-Y (6.8E-08-Y(6.7E-13) ) 


5 


DERIV 


5 


Y=YO+AD*VY 

AP1=    AP2 

ANl  =    AN2 

VX»   VXO+AD*ANl 

VY=    VYO+AD*APl 


RHO 


=2.4E-4-Y(6.8E-8-Y(6.7E-13) ) 


DERIV 


*A  IS  A   IN  PLM 
**  T  IS  TM  IN  PLM 


X 


T=  T+D 

X=  X+D*(VXO+AA*(VX-VXO) ) 
Y=  YO+D* (VYO+AA* (VY-VYO) } 
VX=  VXO+D*(ANl+AA*(AN2-ANl) ) 
VY=  VY0+D*(AP1+AA*(AP2-AP1) ) 
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DERIV   SUBROUTINE 


c 


ENTRY 


) 


GLOBAL 
VARIANCES 


£ 


IR£G=  1 


5. 


V=  SQRT(VX*VX+VY*VY) 
CM=  V* (8.9E-04+3 . 2E-09*Y)+DM 


^Jl 


IREG=  2 


^xr^ 


CKDG=  DKG+CF* (CC ( I REG , 1 , MSTG)  +  (CC(IREG, 2 , MSTG) +CC ( I REG , 3,MSTG)  *CM)  *CM) 


V 


HH=  TH/V-RHO*CKDG*V 
AN2=  HH*VX 
AP2=  HH*VY-G 


SL 


RETURN 
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APPENDIX  B:   PROGRAM  DEFINITION  OF  VARIABLES 


SYMBOL 

PL/M  if 
FORTRAN   varied 


DEFINITION 


A 

AA 

ALT 

AN1,AN2 

AP1,AP2 

CC 

CF 

CFORM1 

CFORM2 

CKDG 

CM 

CT 

DEG 

DEL 

DKG 

DKG1 

DKG2 

DM 

DMAX 

DM1 

DM2 

DS 


Al 


Runge  Kutta  parameter 

0.5/Runge  Kutta  parameter 

Weapon  release  altitude 

Runge  Kutta  variables  (x  direction) 

Runge  Kutta  variables  (y  direction) 

Matrix  of  drag  curve  coefficients 

CKDG  stretch  factor 

CKDG  stretch  factor  for  the  first  stage 

CKDG  stretch  factor  for  the  second  stage 

Bomb  coefficient  times  drag  coefficient 

Mach  number 

Matrix  of  Mach  curve  coefficients 

Weapon  release  angle 

ArcTan  (ejection  velocity/weapon  velocity) 

Shift  in  CKDG 

Shift  in  CKDG  for  the  first  stage 

Shift  in  CKDG  for  the  second  stage 

Shift  in  Mach  number 

Largest  integration  step  size 

Shift  in  first  stage  Mach  number 

Shift  in  second  stage  Mach  number 

Integration  step  size 
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SYMBOL 

PL/M  if 
FORTRAN   varied 


DEFINITION 


DTI 

FN 

FRACT 

G 
HH 

IANG 
I  BOTH 

IDNO 

IREF 

IREG 

IV 

IY 

MSTG 

RAD 

RHO 

SET 

SL 

T 
TH 

THETA 

U 

V 

vx 


TM 


Initial  integration  step  size 

Equal  to  the  Thrust  for  the  first  stage 

Parameter  for  Runte  Kutta  integration 
(value  between  zero  and  one) 

Acceleration  due  to  gravity 

Total  drag  function 

Number  of  release  angles  integrated 

Key  determining  if  two  drag  curves 
are  required  for  the  weapon 

Identification  number  of  each  weapon 

Reference  for  Mach  curve  coefficients 

Mach  region  index 

Number  of  velocities  integrated 

Number  of  altitudes  integrated 

Weapon  stage  index 

Conversion  factor  for  radians 

Air  density 

Key  to  initializing  variables 

Slope  coefficient  for  retarded  snakeye 
fins  deployment  time 

Elapsed  time  of  flight  from  weapon  release 

Thrust/Mass  of  rocket 

Adjusted  weapon  release  angle 

Aircraft  speed 

Velocity  of  weapon 

Velocity  component  of  V  in  the  x 
direction 
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SYMBOL 

PL/M  if 
FORTRAN   varied DEFINITION 

VXA  VX  at  weapon  release 

VXO  VX  at  start  of  present  integration  step 

VY  Velocity  component  of  V  in  the  y 

direction 

VYA  VY  at  weapon  release 

VYO  VY  at  start  of  present  integration  step 

VZ  Muzzle  velocity  of  guns 

X  Weapon  ground  range  from  release 

Y  Weapon  altitude  above  sea  level 

YO  Y  at  start  of  present  integration  step 

YT  Target  altitude  above  sea  level 
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FORTRAN  AND  PL/M  RESULTS 


IDNO 

VKTS 

ALT 

DEG 

FOR^ 

PRAN 

PL/M 

TOF 

RANGE 

TOF 

RANGE 

1 

450. 

3000. 

-10. 

10.26 

7515. 

10.27 

7518. 

2 

450. 

3000. 

-10. 

10.41 

7393. 

10.41 

7395. 

3 

450. 

3000. 

-10. 

10.32 

7467. 

10.32 

7470. 

4 

450. 

3000. 

-10. 

10.73 

7189. 

10.74 

7191. 

5 

450. 

3000. 

-10. 

10.91 

7085. 

10.8  7 

7078. 

6 

450. 

3000. 

-10. 

11.  13 

695  5. 

11.  14 

6957. 

7 

450. 

3000. 

-10. 

10.32 

7466. 

10.33 

7469. 

8 

450. 

3000. 

-10. 

10.59 

7292. 

10.6  0 

7294. 

9 

450. 

3000. 

-10. 

10.29 

7492. 

10.29 

7495. 

10 

450. 

5000. 

-10. 

10.25 

7521. 

10.  26 

7524. 

11 

450. 

3000. 

-10. 

10.25 

7529. 

10.26 

7532. 

12 

450. 

3000. 

-10. 

10.24 

7536. 

10.24 

7  5  39. 

13 

450. 

3000. 

-10. 

10.22 

7547. 

10.23 

7550. 

14 

450. 

3000. 

-10. 

10.31 

7478. 

10.31 

7480. 

15 

450. 

5000. 

-10. 

10.40 

7424. 

10.41 

742  7. 

16 

450. 

3000. 

-10. 

10.  26 

7  516. 

10.27 

7  519. 

17 

450. 

3000. 

-10. 

10.50 

7368. 

10.51 

7371. 

18 

450. 

3000. 

-10. 

16.55 

3852. 

16.57 

3847. 

19 

450. 

300  0. 

-10. 

20.84 

4620. 

20.87 

4619. 

20 

450. 

3000. 

-10. 

14.92 

5771. 

14.94 

5768. 

21 

450. 

3000. 

-10. 

16.83 

5484. 

16.87 

5  4  81. 

22 

450. 

3000. 

-10. 

17.96 

3299. 

17.98 

3296. 

23 

4  5  0. 

5000. 

j 

-10. 

10.67 

10  9  4  0. 

10.69 

10935. 

FORTRAN  AND  PL/M  RESULTS  (continued) 


IDNC 

VKTS 

ALT 

DEC 

FORTRAN 

PL/M 

TOF 

RANGE 

TOF 

RANGE 

24 

450. 

3000. 

-10. 

7.82 

12814. 

7.85 

12815. 

25 

4  50. 

3000. 

-10. 

31.70 

1657. 

31.70 

1656. 

26 

450. 

3000. 

-10. 

36.29 

1384. 

36.30 

1382. 

27 

450. 

3000. 

-10. 

41.58 

1213. 

41.58 

1210. 

28 

450. 

3000. 

...  ..  _  .  .  .. 

-10. 

20.56 

2444. 

20.57 

2442. 
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